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Abstract 

We consider the Cacuhy problem for a viscous compressible rotat- 
ing shallow water system with a third-order surface-tension term in- 
volved, derived recently in the modelling of motions for shallow water 
with free surface in a rotating sub-domain [T3] . The global existence of 
the solution in the space of Besov type is shown for initial data close 
to a constant equilibrium state away from the vacuum. Unlike the 
previous analysis about the compressible fluid model without coriolis 
forces, see for instance [§J [T3], the rotating effect causes a coupling 
between two parts of Hodge's decomposition of the velocity vector 
field, and additional regularity is required in order to carry out the 
Friedrichs' regularization and compactness arguments. 

1 Introduction 

The nonlinear shallow water equation is used to model the motion of 
a shallow layer of homogeneous incompressible fluid in a three dimensional 
rotating sub-domain and, in particular, to simulate the vertical average dy- 
namics of the fluid in terms of the horizontal velocity and depth variation. In 
general, it is modeled by the three dimensional incompressible Navier-Stokes- 
Coriolis system in a rotating sub-domain of M 3 together with a (nonlinear) 
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free moving surface boundary condition for which the stress tension is evolved 
at the air-fluid interface from above and the Navier boundary condition of 
wall-law type holds at the bottom. Under a large-scale assumption and 
hydrostatic approximation, the nonlinear shallow water equation has been 
derived recently in [T3], [18] . Usually, the nonlinear shallow water equations 
take the following form of compressible Navier-Stokes equations 

h t + div(toi) = 0, 

(hu) t + div(/m ® u) + ghVh + f(hu) 1 - 

= div(2f (fc)£>(u)) + V(A(/i)divu), !1 " U 
k h(0) = h , u(0) = u , 

where h(t,x) is the height of the fluid surface, u(t, x) = {u l {t,x),u 2 {t,x)) T 
is the horizontal velocity field, u ± (t, x) = (— u 2 (t, x), x)), x = (xi,x 2 ) G 
1R 2 , D(u) = |(Vu + (Vu) T ), g > is the gravitational acceleration, f > is 
the Coriolis frequency, £ ^ and A are the dynamical viscosities satisfying 

For the shallow water system (11. ID . there is a mount of work to deal with 
the global well-posedness of strong solutions subject to some small initial 
perturbation of a constant state or the global existence of weak solutions 
for large initial data. When the viscosities satisfy £(h) = h and A = 0, 
and the effect of the Coriolis force and/or third-order surface tension term is 
omitted (/ = 0, (3 = 0), the local existence and uniqueness of classical solu- 
tions to the Cauchy-Dirichlet problem for the shallow water equations with 
initial data in C 2+a was studied in [I] using Lagrangian coordinates and 
Holder space estimates. Kloeden and Sundbye [161 E2] proved the global ex- 
istence and uniqueness of classical solutions to the Cauchy-Dirichlet problem 
using Sobolev space estimates by following the energy method of Matsumura- 
Nishida [19]. Sundbye [23] proved also the existence and uniqueness of clas- 
sical solutions to the Cauchy problem using the method of [19]. Wang-Xu, in 
|24j . obtained local solutions for any initial data and global solutions for small 
initial data h — h , u G H 2+S (M?) with s > 0. The result was improved by 
Chen-Miao- Zhang and Haspot to get global existence in time for small initial 
data ho — ho G S° 1 fl B\ 1 and u G x as a special case in [7J [15]. Cheng- 
Tadmor discussed the long time existence of approximate periodic solutions 
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for the rapidly rotating shallow water for initial data (h ,u ) e H m (T 2 ) 
with m > 5 where the viscous terms are absent (i.e. £ = A = ) in [S]. 
The global existence of weak solutions for arbitrarily large initial data is 
established in one dimension [17], where the vanishing of vacuum states in 
finite time is shown, and in multi-dimensional bounded domain with spher- 
ical symmetry [14J with the help of the Bresch-Desjardins entropy [2] and 
the L -stability compactness argument [20]. The global existence of weak 
solutions for arbitrarily large initial data is shown by Bresch-Desjardins [3] 
where additional drag friction and capillary terms are involved to construct 
an global approximate solutions. An related systems with a third-order term 
stemming from the capillary tensor also have been considered by Danchin- 
Desjardins for a compressible fluid model of Korteweg type p2] with constant 
viscosity coefficients, and the global existence of strong solution is shown. 

In the present paper, we consider the global existence of the Cauchy 
problem for the 2D viscous shallow water equations 

' ht + div(/m) = 0, 

(hu) t + div(/m <g> u) + ghVh + f(hu) x 
< (1.2) 
= 2//div(/iD(u)) + 2/iV(/idivu) + (3hVAh, 

k h(0) = h , u(0) = uo, 

which corresponds to ( 11. ip for the case 2£(h) = X(h) = 2fih with /i>0a con- 
stant, and is derived recently in [18] with a third-order surface tension term 
involved by considering second order approximation and parabolic correction 
where (3 > is the capillary coefficient. Although there are many mathe- 
matical results about the shallow water equations (11.11) . there is no analysis 
about Eq. ( 11. 2ft . It also should be mentioned that the global existence of 
weak solutions does not apply here since the Bresch-Desjardins entropy [2] 
is not satisfied for Eq. (11.21) . In addition, the classical theory does not cover 
the case with coriolis force and capillarity term involved. 

We investigate the global existence of strong solution in some Besov space. 
Although we also make use of the Hodge's decomposition to separate the ve- 
locity field into a compressible part and an incompressible part, unlike p,[TU] 
we finally obtain a coupled system due to the rotating effect of the Corio- 
lis force. In fact, it can not be decoupled into a system involving only the 
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compressible part and a heat equation containing only the incompressible 
part because of the appearance of the Coriolis frequency, which leads to a 
strong coupling between the gradient vector field part and divergence free 
part of the fluid velocity in terms of the Hodge's decomposition. Thus, we 
have to investigate the whole system of the height, the compressible velocity 
field part and the incompressible velocity field part. With the help of the 
Littlewood-Paley analysis and hybrid Besov spaces, we obtain the a priori 
estimates in Chemin-Lerner type time-spatial spaces which are necessary in 
order to use the interpolation theory of time-spatial spaces involving hybrid 
Besov spaces. Then we use a classical Friedrichs' regularization method to 
construct approximate solutions and prove the existence of a solution by 
compactness arguments. For the uniqueness of solutions, due to the contri- 
bution of the third-order surface tension term, we can prove it in a larger 
space than that for the existence and we do not need more regularity on the 
spaces. 

For the convenience of the statement of main results, we note that B^f 2 
is a hybrid Besov space defined in the next section, the space E s is defined 
by 

E s ={(/», u) G C ([0, oo); B^) H L 1 (o, oo; B*+ 3 ' s+2 ) 

x (C ([0, oo); S'i 1 ) n L 1 (0, oo; i^ 1 ))' }, 

and C([0, oo); B^f 2 ) is the subset of functions of the Chemin-Lerner type 
space L^^B^Y' 2 ) defined in the next section which are continuous on [0, oo) 
with values in Bf^f 2 ■ 

For the initial data h , we suppose that it is a small perturbation of some 
positive constant h . The main theorem of this paper reads as follows. 

Theorem 1.1 Let e G (0, 1), h - h G £° i +e and u G B^'J. Then, there 
exist two positive constants a small enough and M such that if 

\\ho ~~ ho\\ 60,i+e + ||uo|| rO,e ^ a, 

then (II. 2p yields a unique global solution (h,u) in (h , 0) + (E 1 nE 1+e ) which 
satisfies: 

\\(h - ^ Q ,u)|| E i nE i+ e ^ M(\\h - h \\ S o,i+e + Huoll^e), 
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where M is independent of the initial data. 

The paper is organized as follows. We recall some Littlewood-Paley the- 
ories for homogeneous Besov spaces and give the definitions and some prop- 
erties of hybrid Besov spaces and Chemin-Lerner type spaces in the second 
section. In Section 3, we are dedicated into proving of the a priori estimates. 
In Section 4, we prove the global existence and uniqueness of solution for 
small initial data by using a classical Friedrichs' regularization method and 
compactness arguments. 

2 Littlewood-Paley theory and Besov spaces 

Let ip : R 2 — > [0, 1] be a radial smooth cut-off function valued in [0, 1] 



Thus, ip is supported in the ball {£gR 2 :|£|^4/3}, and (p is also a smooth 
cut-off function valued in [0, 1] and supported in the annulus {£ : 3/4 < |£| ^ 
8/3}. By construction, we have 



such that 



1, If I ^ 3/4, 

smooth, 3/4 < |f | < 4/3, 
0, |f | ^ 4/3. 



Let </?(£) be the function 



^(0:=^/2)-V(0- 



5>(2" fc O = l, V^O. 



One can define the dyadic blocks as follows. For k G Z, let 



A fc / := ^-V(2- fe 0^/- 



The formal decomposition 




(2.1) 
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is called homogeneous Littlewood-Paley decomposition. Actually, this de- 
composition works for just about any locally integrable function which yields 
some decay at infinity, and one usually has all the convergence properties of 
the summation that one needs. Thus, the r.h.s. of ( 12. ip does not necessarily 
converge in ^'(R 2 ). Even if it does, the equality is not always true in ^'(R 2 ). 
For instance, if / = 1, then all the projections A^f vanish. Nevertheless, 
( 12.11) is true modulo polynomials, in other words (cf. [TH [21]), if / G ^'(R 2 ), 
then Efcez A fc/ converges modulo ^[R 2 ] and fl2ZD) holds in ^'(R 2 )/^[R 2 ]. 

Definition 2.1 Let s G R ; 1 < p, q ^ oo. For f G y'(R 2 ), we write 

ii/ii^ = E 2fa ii A */i^- 

fcez 

A difficulty comes from the choice of homogeneous spaces at this point. 
Indeed, || • cannot be a norm on {/ G o$^'(R 2 ) : H/H^ < °°} because 
II/IIbIj = means that / is a polynomial. This enforces us to adopt the 
following definition for homogeneous Besov spaces (cf. [TO]). 

Definition 2.2 Let s G R and m = —[2 — s]. If m < 0, then we define 
^(R 2 ) as 

Ki = {/ e ^'( R2 ) : Wf\\B h < oo and f = J2 A *f m ^'(K 2 )}. 

fcez 

If m ^ 0, we denote by ^ m the set of two variables polynomials of degree 
less than or equal to m and define 

Ki = {f e y'^l&m : \\f\\ Bil <^andf = ]T A k f in ^'(R 2 )/^ m }. 

fcez 

For the composition of functions, we have the following estimates. 

Lemma 2.3 (pH Lemma 2.7]) Let s > and u G B s 2 l n L°° . 

i) Let F G W£J +2 '°°(R 2 ) such that F(0) = 0. Then F(u) G B s 21 . More- 
over, there exists a function of one variable Cq depending only on s and F, 
and such that 

\\F(u)\\ H ^C (\\u\\ L ^)\\u\\ Bii . 
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ii) Ifu,v E B\ x , (v - u) G B s 21 for a s G (-1, 1] and G G W^(R 2 ) 
satisfies G'(0) = 0, then G{v) — G{u) G B 2l and there exists a function of 
two variables C depending only on s and G, and such that 

\\G(v) - G(u)\\zs 2 i s= C7(||u||r«, (iMIe^ + IMI^J Ik ~ "11^- 

We also need hybrid Besov spaces for which regularity assumptions are 
different in low frequencies and high frequencies [TU]. We are going to recall 
the definition of these new spaces and some of their main properties. 

Definition 2.4 Let s, t G R. We define 

ll/b- = £ 2 *1 A */IU a + J> fc lA fc /|| L2 . 

Lei m = — [2 — s], we i/ien define 

5*j(R 2 ) = {/ G J^'(M 2 ) : ll/H^j < oo} , i/m < 0, 
^(R 2 ) = {/ e y'(R 2 )/^ m : ||/||^,t < oo} , ifm > 0. 

Lemma 2.5 VFe /iai>e i/ie following inclusions. 

(i) We have B 2 'l = B S 2V 

(ii) Ifs^t then B 2 j = B s 21 n B l 2l . Otherwise, B 2 '^ = B s 2l + B\ x . 

(Hi) The space B 2 'l coincides with the usual inhomogeneous Besov space 

B 2,l- 

(iv) lfs x ^ s 2 and tt ^ t 2 , then B^f 1 B s 2 f 2 . 

Let us now recall some useful estimates for the product in hybrid Besov 
spaces. 

Lemma 2.6 ( [TOl Proposition 2.10]) Let s h s 2 > and f,g G L°° n 

S^ S2 . T/ien fg G L^f 2 and 

ll/^lls^ < ||/||L°°lblls 2 y2 + 11/11^3 ||^|U«. 

Lei si, s 2 , ii, *2 ^ 1 snc ^ mm ( s i + s 2,ti + £2) > 0, / G L?^ 1 and 
g G L?0 t2 . T/ien fg G S^-iA+^-l flnd 

||/fi'||^i+ s 2-i,*i+'2-i < ||/HBn.*i||o||^2,t2. 
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In the context of this paper, we also need to use the interpolation spaces of 
hybrid Besov spaces together with a time space such as L p (0, T; B^'i)- Thus, 
we have to introduce the Chemin-Lerner type space which is a refinement of 
the space L p (0, T; B^'i)- 

Definition 2.7 Let p G [1, oo], T G (0, oo] and s±, s 2 G M. Then we define 

Wf\\l p T {B s 2 'X) = S ^ 2kS \\^kf\\LP{0 ) T;Ll) + ^ Afc/ 1| Lr(0,T;L 2 ) ■ 

k^O k>0 

Noting that Minkowski's inequality yields ||/HiP(jj|.*) ^ ll/llz£,(i?2'*)> we 
define spaces L^B^'i) as follows 

L P T (B s 2 i) = {/ G L P T (B%) : Wfh^ < oo} . 

If T = oo, then we omit the subscript T from the notation L^B^'i), that 
is, LP(B2i) for simplicity We will denote by C([0, T]; B^'i) the subset of 
functions of (B21) which are continuous on [0, T] with values in 

Let us observe that L^B^i) = L^B^), but the embedding L^B^'i) C 
L^B^'i) 1S strict if p > 1. 

We will use the following interpolation property which can be verified 
easily (cf. pQ). 

Lemma 2.8 Let s, t, si, t±, S2, ti G M and p,pi,p2 6 [1, oo]. We have 

where ± = £ + s = flsj + (1 - #)s 2 and t = ^ + (1 - 6)t 2 . 

Now, we define the following work space. 
Definition 2.9 For T > and s 6 1, we denote 

E^j^ulGC^T^^jni 1 (0,T;5 2 S + 3 < S+2 ) 

x (C ([0, I*]; J3£ l ) H (0, T; 2 } 

and 

\\(h,u)\\ E *. =\\h\\ i¥{§ s-i, S) + Hull^^-i) 

+ INIz4(^+ 3 - 3 + 2 ) + II u IIl1(b5+ 1 )- 
We use the notation E s if T = +oo, changing [0, T] into [0, +oo) in the 
definition above. 
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3 A priori estimates 



Noticing that divZ)(u) = |Vdivu + |Au and substituting h by h + ho 
in (II. 2D . we have 

' h t + u • Vh + ho divu = — /i divu, 
u ( + U'Vu - fiAu - 3/iV divu + ftr 1 + gV/i - 13V Ah 

VhD(u) + V/idivu (3.1) 
h + fiQ 

k /i(0) = /i - h , u(0) = u . 

For all s G 1, we denote A 5 / = ^" _1 (ICr/)- Le t c = A -1 divu and 
d = A -1 div^u where div^u = V -1 • u and = (— <?2, c\). Then, it is easy 
to check that 



u 



-A _1 Vc - A _1 V- L d. 



Now, we can rewrite the system (13.11) in terms of these notations as the 
following: 



where 



/it + u • Vh + h Ac = F, 
q + u ■ Vc - 4/xAc -id- gAh - pA 3 h = G, 
dt - //Ad + fc = A" 1 div^H, 
u = -A _1 Vc- A" 1 V ± rf, 
/i(0) = /i - ^o, u(0) = Uo, 

F = — /idivu, 

G = u • Vc + A" 1 divF, 

„ VMXu) + V/idivu 
H = -u ■ Vu + 2// V — 7 • 

h + h 

For these equations, we study the following system: 

h t + v ■ Vh + /i Ac = F, 
ct + v • Vc - 4/xAc - fd - gA/i - /3A 3 /i = G, 
d t — fiAd + ic = P, 



(3.2) 



(3.3) 
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where v is a vector function and we will precise its regularity in the following 
proposition. 

Proposition 3.1 Let (h,c,d) be a solution of (13.31) on [0,T), T > 0, < 
s ^ 2 and V(t) = J * ||v(r)||_g2 dr. The following estimate holds on [0, T): 

+ / (ll^)ll^+3, s+ 2 + || C (r)||^+i + ||d(r)||^+i)dT 

** 

< Ce™®(\\h{0)\\sg* + 110(0)11^-! + |K0)b ril 

+ f e~ CV{T) (\\F(t)\\ § s-i„ + \\G(r)\\^ + ||P(r)||^-i) dr), 

«/ 

where C depends only on s, ho and coefficients fi, f, g and f3. 
Proof. Let (h, c, d) be a solution of (13. 3p and we set 

(7i, c, J, F, G, P) = e"^ y(t) (h,c,d,F, G, P) . 
Thus, ( 13. 3 p can be transformed into 

ht + v ■ Vh + h A5 — F — KV'(t)h, 

c t + v • Vc - 4/iAS - fd - gAh - (3A 3 h = G- KV'(t)c, (3.4) 
d t - fi Ad + fc = P - KV'(t)d. 

Applying the operator to the system (13.4p . we obtain the following sys- 
tem by noting (h k , c k , d k , F k , G k , P k ) = (A k h, A k c, A k d, A k F, A k G, A k P): 

d t h k + A fc (v • Vh) + h A5 k = F k - KV'{t)h k , 

d t c k + A fc (v • Vc) - 4/iAgfc - fd fc - gA/i fc - (3A 3 h k = G k - KV'(t)c k , 
8 t d k - nAd k + ic k = P k - KV\t)d k . 

(3.5) 

To begin with, we consider the case where v = 0, K = and F = G = 
P = which implies that (13. 5p takes the form 

d t hk + hoAc k = 0, 

d t ~c k - 4/iAc, - f4 - gA/i fc - /?A 3 /i fc = 0, (3.6) 
d t d k - /iA4 + fc fe = 0. 
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3.1 The case of high frequencies 



Taking the L 2 scalar product of the first equation of (13.61) with h k , of the 
second equation with c k , and the third one with d k , we get the following two 
identities: 

' 1 d „7 



2dt 
1 d 

2di 



\hk\\ L 2 + h (Ac k , h k ) = 0, 

\ck\\ 2 L 2 +4//||Ac fe |||2 -f(d k ,c k ) -g(Ah k ,Ck) 
-P(A 2 ~h k ,Ad k )=0, 



(3.7) 



^IMfcllia + n\\Ad k \\ 2 L 2 +f(c fc ,4 



0. 



Now we want to get an equality involving Ah k . To achieve it, we take L 2 
scalar product of the first equation of (13.61) with A 2 h k and Ac k respectively, 
then take the L 2 scalar product of the second equation with Ah k and sum 
with both last two equalities. This yields 

1 d 



— r i \Ah k \\U 
d 



+ h (A c fc , Ah k 



0. 



— (Ah k , c fc ) +/io||Acfc||i2 -i(dk,Ah k ) -g\\Ah k \\ L 2 
- (3\\A 2 hk\\ 2 L 2 + M^c k ,A 2 ~hk) = 0. 
Let K\ > be a constant to be chosen later and denote for k > 

a l = T-\\hk\\ 2 L2 + "T~II a Ml 2 + llcfellia + \\dk\\ 2 L 2 - 2Ki(Ah k , c k ) 
h h 

By a linear combination of (13.71) and (13. 8ft . we can get 

1 d 



(3.1 



2dt 



at + (4/i - /lo^OIIAcfcHi, + g^HA^Hi, + pK^A'hkWt* 
+ n\\Ad k \\ 2 L z - AfiKi(Ac k , A 2 h k ) +iK 1 (d k ,Ah k ) = 0. 



(3.9) 



Using Schwartz' inequality, Young's inequality and Bernstein's inequality 

4 



\Ad k \\L 2 , 

we find, for any positive numbers Mi, M 2 , M 3 , that 

1 



(Aq..A-// a ) ^||Ar A ||i, + — ||A-// "- 



2Mi 



fcllL 2 ) 



12 



C.C.Hao, L.Hsiao, H.-L.Li 



(d k ,Ah k 



1 



9 



k\\L 2 



2M 2 

1 



W^hkWh, 



(Ah k ,c k ) <^ — \\Ah k \\ 2 L2 + 2]f 



2 

L- ' 



Thus, we need to determine the values of Ki, M±, M 2 and M 3 such that 

*-^"Wf >0 ' /J -2^ >0 ' g -2lS >0 ' 

/i ^>0, f-KiM 3 >0, l-v^>0. 

9 h M 3 

One can verify that the above inequalities will hold if one has 



< K\ < min 



4/i/? 2 /3 9/ig 



/i /3 + 5/i 2 ' 3 V V 4f 2 



(3.10) 



5u _ f 9u w 2 
M x = — , M 2 = — + — — , M 3 = -W=-. 
2/3' 4g 16fKi' 3 3y/i 

Hence, we obtain 

Therefore, there exists a constant c > such that 

i|<4 + «j»o4<o. 

In the general case where F, G, P, .K" and v are not zero, we have, with 
the help of Lemma 6.2 in pU], that 



a 



I + (c2 2k + KV) 



1 d 
2~dt 

g ~ ~ 

fc fc ) + (G fc , c fc ) + (P fc , c fc ) 

+ £-(AF k ,A~h k ) - Ki(AF k , c fc ) - K^Ah) 

- f (A fc (v ■ Vfc), - (A fc (v ■ Vc), c fc ) - |-(AAfc(v • V/l), A/i fc ) 
fto «0 

+ K x (AAfc(v • Vfc), c fe ) + Ki(A fc (v • Vc), A/ifc) 
<Xk(\\Fk\\v + ll^fclk 2 + 11-PfcllLJ" + IIAFfclb 



< 
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+ 7fc2" 
+ 7fc2 



-fe(s-l) I 



-k(s-l) I 



R2 



V R2 



ds-1 
^2,1 



7*2-*^ 
3^+7*1^11^(2 



V £)2 C nS-1 

-*(«-!) ||S|| . . 
II ll_B 2 



_l_ 2~fc(s-i) 



where ^ fc 7fc ^ 1 and s G (0,2]. 



3.2 The case of low frequencies 

We replace the second equation of (13.81) by the following equation 



d_ 

It 



(A 3 h k , c k ) + h \\K 2 ~c k \\ 2 L 2 - f(A4, A 2 h k ) - g\\A 2 h 



k\\L 2 



(3.11) 



- P\\A 6 h k \\ 2 L2 + Afj,(A 2 c k ,A 6 h k ) ^ 



Let Ko > be a constant to be chosen later and denote for k ^ 



«2 = f \\h\\h + r\\M h \\h + 

h h 



7, 1 1 2 



A linear combination of (13.71) . the first equation of 
1 d 



and (13. lip yields 



~o2 + 4//||A5 fc ||i 2 - ^ 2 ||A 2 5 fc ||i 2 + gK 2 \\A 2 h k \\ 2 L2 



(3.12) 



+ /?ir 2 ||A 3 /i fc ||| 2 + /i||A4||i2 - 4/xK 2 (A 2 4, A 3 /i fc ) 
+ fK 2 (Ad k ,A 2 h k ) = 0. 

Using Schwartz' inequality, Young's inequality and Bernstein's inequality 

||A/i fe || L 2 <C ^2 k fh k \\ L2 , 
we find, for any positive numbers M 4 , M5, Me, that 



(A 2 g fc ,A^ fc ; 



8 2 M 4 „ ,. „, I 
2-3 2 



l|Ac fe ||i 2 + — — ||A^ 



2M 4 



2 



(Ad k ,A 2 h k ) ^\\Ad k \\ 2 L2 + ^H|A 2 Mi 2 , 



(A 3 /i fc ,g fc ) 



<- 



8 4 M fi 



\\Ah k \\ 2 L 2 + 



I ~ l|2 
|Cfe|| L 2. 



2.3 4 "'"' ,|IJ/ ' 2M 6 
Thus, we need to determine the values of K 2 , M 4 , M 5 and M 6 such that 



8% , T ^ 8 2 M 4 

4// - -/i iT 2 - 2732 



4u f 
>°> / ? -7TT^>0, g~— >0, 



2M 4 



2M fi 
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fK 2 M 5 8 8 4 K 2 

u > 0, f- 7 K 2 M 6 > 0, 1 - > 0. 

P 2 ' /i 2-3 4 ' M 6 

One can verify that the above inequalities will hold if one chooses 



< K 2 < min 



/ 3 2 ■ Afi/3 3 2 [J 4/ig 



M 4 = -g, M 5 = — + M 6 = -Wf . 
2p 4g f A 2 8^ y ft 

Hence, we obtain 

csctfe ^ \\hk\\h + l|AM£a + Hcfcll^a + \\dk\\ 2 L 2 < C4C& (3.13) 
Therefore, there exists a constant c > such that 

In the general case where F, G, P, K and v are not zero, we have, with 
the help of Lemma 6.2 in pU], that 



al + (c2 4fe + KV) a 2 k 



1 d 

+ £-(AF k ,A~h k ) - K^F^Ck) - Ki(G k ,A 3 h k ) 

- f (A fc (v ■ Vh), h k ) - (A fc (v ■ Vc), c fc ) - |"(AA A (v ■ V/i), A/i fc ) 

+ iTi(A 3 A fc (v ■ Vfc), c fc ) + #i(A*(v • Vc), A 3 /i fe ) 
<a fc (||F fe || L 2 + ||(5 fc || L 2 + ||P fe || L 2 

+ 7!b2-*^||v||^J|/i||^-i + 7 fc 2-^- 1 )||v||^ ||c||^ 



+ l k 2-^\\v\\^Jh\\^ +7^1^11^(2-^11511^-, 
+ 2-*<- 1 )||/i||^i)), 

where J2k 7fc ^ 1 an d s e (0, 3]. 
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Thus, combining two case of high and low frequencies, we obtain for any 
k G Z 



1 d 



~—a\ + (c2 2k min(l, 2 2k ) + KV') 



Q 



<a k (\\F k \\ L . + \\G k \\ L 2 + \\P k \\ L2 + ||AF fc || L2 (3.14) 
+ 7*2-^^11(^5)11^-1,.^), 

where we choose c = min(c, c) and V(£) = J ||v||_^a . 

We are now going to show that the inequality (13.141) implies a decay for 
h, c and d. 

3.3 The damping effect for h 

Dividing (I3.14p by a k , we get 

ja k (t) + (c2 2k min(l, 2 2k ) + i^V) a, 

<\\h\W + H^Hl* + II^IU 2 + IIAFfcllia + T^-^y'llft^ll^-i,.^-!. 

(3.15) 

Integrating over [0, £], we have 

a fc (t) + c2 2fc min(l,2 2/c ) / a fe (r)dr 
<a fc (0) + f (\\F k (r)\\ L , + ||G,(r)|| L2 + \\P k (r)\\ L , + ||AF fc (r)|| i2 W 



o 



+ / V(t) 







7fe (r)2- fc (- 1 ) || (fc, 5) || ^..^-i - Ka k {r) 



dr. 



By the definition of a 2 , we have 

2*< S - 1 W ^2 fc *max(l,2- fc )||/i fc || L2 + 2 k ^\\~c k \\ L 2 
+ 2 fc ( s - 1 )||4|| i2 , WkeZ. 

Thus, we have, by taking K large enough, that 



(3.16) 



(3.17) 



^ 0. 
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Changing the functions (h, c, d, F, G, P) into the original ones (h, c, d, F, G, P) 
and multiplying both sides of (13.161) by 2 k ( s ~ 1 \ According to the last inequal- 
ity and due to (13.161) and (13.171) . we conclude after summation on k in Z, 
that 

INIl-^- 1 '*) + IMIl-gb^ 1 ) + Nllz-^ 1 ) + c / \\H T )\\B s + 1 s ' s+2dT 

J 

+ cJ2 I 2 fc(s+1) min(l, 2 2fc ) || Cfe (r) W^dr 

k&L ^° 

+ cJ2 / 2 fc ( s+1 Vin(l,2 2fc )||4(r)|| L 2rfr ( 3 ' 18 ) 
<e^W||(MO),c(0),rf(0))b 2S -M x( ^- 1)2 

+ e cv W £ e-^W || (F, G, P)(r) || ^-m^-^t. 

3.4 The smoothing effects of c and <i 

Once the damping effect for h is established, it is easy to get the smooth- 
ing effect on c and d. Since (13.181) implies the desired estimate for high 
frequencies, it suffices to prove it for low frequencies only. We therefore 
suppose in this part that k ^ 0. 

Taking the L 2 scalar product of the last two equations of f !3. 51) with c k 
and dk respectively, we have 

' 1 d ~ 

oTJillMla + h (c k ,Ah k ) 



2dr 



(F k , h k ) - KV'{t)\\h k f L2 - (A fe (v • Vh), h 



k) 



\j t \\~ c k\\h + 4/i||Acfc|||2 - f(4, ~c k ) - g(A4, Ck) - (3(A 3 h k , 4) (3-19) 

= (G k ,c k ) - KV\t)\\c k \\v - (A fc (v ■ Vc),c fc ), 
\j t \\dk\\h + fJ>\\M k \\ 2 L 2 + f(4,c fe ) = (Pfc,4) - W(t)||4||i 2 . 

Define 6 12 , = ^||4 11^2 + 11^11^2 + ||4|li2- By using Lemma 6.2 in [TP] , 
(I3.19P yields, for a constant c > 0, that 

^ 2 + C2 2k 9 2 k <^.(||A 3 ^|| L 2 + \\G k \\# + \\P k \\v>) 
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+ 6 k V'(t)(C lk 2- k< - s - 1 \\\5\\ B s- 1 + \\h\\ H -i) -K9 k ). 



Dividing by 9 k and integrating over [0,t], we infer 



t rt 
2k i 



9 k (t)+C 2 2k 9 k (r)dr^9 k (0) + C [\\G k {r)\\ L 2 + \\P k {r)\\ L ,]dr 
Jo Jo 

+ C /V'(r) Tfe (r)2-^- 1 )(||c(r)|| 4s - 1 + fh(r)\\^)dr. 



o 



Therefore, changing the functions (h, c, d, F, G, P) into the original ones, we 
get 

^2^- 1 )||/. fc (t)|| z§?{L2) + ^2^- 1 )||c fc (t)|| z§?(L2) 



T ' 

fc<0 fc<0 



k^O ^° k^O 

+ C f Y,^ {s+ " ) \\dk{r)\\ L ,dr 
<e cv ^ || (MO), c(0), d(0)) || ^-m x( ^- 1)2 + f e cv ^ || (G(r), P(r)) || ^dr 

J 

+ e cvit) (\\c\\ Lr{Kil) + \\h(T)\\ i¥{§ s-i, S) ) 
Using (13.181) . we eventually conclude that 

C fY, 2k{S+l) \\ c ^)\\vdT + C r^2^ +1 )||4(r)|| L2 rfr 

<e c ^(||(M0),c(0),40))b ri , Sx( ^ ril)2 

+ £ e~ CV{T) || (F(r), G(r), P(r)) || ^ x( ^-i )8 dr) 
Combining the last inequality with (13.181) . we complete the proof. □ 

4 Existence and uniqueness 

This section is devoted to the proof of the Theorem 11.11 The principle 
of the proof is a very classical one. We shall use the classical Friedrichs' 
regularization method, which was used in [5J El El EES] for examples, to con- 
struct the approximate solutions (/i n ,u n ) ne N to (13.11) . and then we will use 
Proposition 13. II to get some uniform bounds on (h n , u n ) n€ N. 
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4.1 Construction of the approximate sequence 

To this end, let us define the sequence of operators (J n )neN by 

and consider the following approximate system: 



h n t + J n (J n u n ■ VJ n h n ) + h AJ n c n = F n , 

c n t + J n ( J n u n ■ VJ„c n ) - A^AJ n c n - f J n d n - gk.J n h n - (3A 3 J n h n = G n , 
<T t - fiAJ n d n + iJ n c n = JnA' 1 div ± H n , 
u n = -A~ x Vc n - A-V^cT, 

(h n ,C n ,d n )(0) = (^A-MivUn^-MivV), 

(4.1) 

where 

h n = J n (ho — ho), u„ = J n Uo, 

F n = -J n (J n h n dwJ n u n ), 

G" = J n {J n u n ■ VJ n c n ) + JnA- 1 divH n , 

ii = -J n u ■ VJ„u + 2/i 77TT^m~\ ' 

C( J n h n + h ) 

with ( a smooth function satisfying 



{ h /A, \s\ ^ ho/ '4, 

s, V2 |s| ^ Zho/2, 

7ho/A, \s\ > 7^ /4. 

smooth, otherwise. 



We want to show that (14.11) is only an ordinary differential equation in 
L 2 x L 2 x L 2 . We can observe easily that all the source term in (14.1 j) turn 
out to be continuous in L 2 x L 2 x L 2 . For example, we consider the term 
JnA' 1 div v ^j" fe ^"^"" . By Plancherel's theorem, Hausdorff- Young's inequal- 
ity and Holder's inequality, we have 

II 7 A -1A- VJ„/l"divJ„u" ! VJ n /i n divJ n u" 
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V T h n div T u n 

£11 : : — L2 £ ||VJ„/i n dlvJ„u"|| L 2| 



C( J„/i n + ^ 







((J n h n + h ) 



\L°° 



4 An 
^ r ||VJ n /i n || L <»||divJ n u n || £2 £ |e|l B( i n) ^ n |Ui||u n || L2 
ho ho <n ' ' 

4r? 3 

£-^ll" IU 2 II U IU 2 - 

Thus, the usual Cauchy-Lipschitz theorem implies the existence of a strictly 
positive maximal time T n such that a unique solution exists which is contin- 
uous in time with value in L 2 x L 2 x L 2 . However, as J 2 = J n , we claim that 
J n (h n ,c n ,d n ) is also a solution, so uniqueness implies that J n (h n ,c n ,d n ) = 
(h n , c n , d n ). So (h n , c n , d n ) is also a solution of the following system: 

h n t + J n {u n ■ Vh n ) + h Q Ac n = F?, 

c n t + J n (u n • Vc n ) - 4/iAc n - fcT - gAh n - (3A 3 h n = G™, 
d n t - fiAd n + fc n = A" 1 div x H?, 
u n = -A _1 Vc n - A^V^d", 
(/*W")(0) = (/^A^divu^A^divV), 



(4.2) 



with 



— Jn{ho — h ), u„ — J n u , 
i? = -J„(/i n divu n ), 

= J n (u n ■ Vc n ) + J^ 1 divfl?, 

V/i"£>(u n ) + V/i n divu r ' 



-u n • Vu r 



2/x- 



C(/i" + /i ) 

The system (14.21) appears to be an ordinary differential equation in the space 

a e L 2 (R 2 ) : supp^a C B(-,n) 

n 

Due to the Cauchy-Lipschitz theorem again, a unique maximal solution exists 
on an interval [0, T*) which is continuous in time with value in L 2 n xL 2 n x L 2 n . 



4.2 Uniform bounds 

In this part, we prove uniform estimates independent of T < T* in E\ D 
E^ +e for (h n ,u n ). We shall show that T* = +oo by the Cauchy-Lipschitz 
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theorem. Denote 

E(0) :=\\ho - hoWtyi+e + ||u ||^o.,, 
E{h,u,t):=\\{h,u)\\ El + \\{h,u)\\ El+ e, 

f n :=sup{t G [0,T*) : E(h n ,u n ,t) ^ ACE(0)} , 

where C corresponds to the constant in Proposition 13.11 and A > max(2, C~ l ) 
is a constant. Thus, by the continuity we have T n > 0. 

We are going to prove that T n = T* for all n e N and we will conclude 
that T* = +oo for any nGff. 

According to the Proposition 13.11 and the definition of (h n ,u n ), the fol- 
lowing inequality holds 

\\(h n ,u n )\\ E , ^""""^^(ll^-SoH^j + 

+ ii^rii^^) + ii u " • v ° n 11^(4%) + n^riLi^))- 

Therefore, it is only a matter to prove appropriate estimates for F™, 
and u n ■ Vc™. The estimate of F™ is straightforward. From Lemma [2.61 we 
have 

W F i\\ L}r{ Bli) < C|l^ll^(^)l|u n L M 4 V) ^ CE\h n ,u n ,T). (4.3) 
With the help of Lemma 12.61 and interpolation arguments, we have 



^llul^ojlul^^) (4.4) 
^CE 2 (h n ,u n ,T). 

In the same way, we can get 

\\n n ■Vu n \\ LUB o i) <:CE\h n ,u n ,T). (4.5) 
To estimate other terms of H±, we make the following assumption on E(0): 

2dACE(0) ^ ho, 
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where C\ is the continuity modulus of B\ x C L°°. If T < T n , it implies 
II /ill l- ^ Cill^lls^ < CiII^IIb^ < CiACE(0) <: -ho. 



Thus, we have 



which yields 



\\h n \\L°°([0,T]xR 2 ) ^ 2^°' 



h n + h E [h , ho] and Q{h n + h ) = h n + h . 
From Lemma 12.61 and 12.31 and interpolation arguments, we have 

vh n ■ Vu n „ 



h 



h n Vh n ■ Vu" 



^\\Vh n ■ Vu^l^^o^ + II - hQ + h n 11^(4%) 

^C , ||V/i n || L oc ( Bo i) ||Vu"||^ ( ^ i) + C\\ =- - ^ 11^(^)11 Vu n || L i ^ 

^cII^IIl-^dII^IIl^b^)^ + ll^ o + ^„lli§?(^, 1 ) / 

^ C7 ll /in |lL ? ?(B 2 ;i 1 )ll Un |lL^(B| il )( 1 + ll^llx??^), 

<:CE\h n , u n , T)(l + £(/i n , u n , T)). 



Similarly, we can get 
XJh n ■ D(u n \ 

ho\\ - h + hn LubIJ < ^ 2 (^ ; u",T)(l + ff(fe",u",T)). (4.7) 
Hence, from fl4.4p - fl4.6l) . we gather 

K ■ Vc n \\ l i t( bo a) + H^l Hli^) 
<C(1 + 4^o x (l + £(/i n , u n , T)))E 2 (h n , u", T). 

Similarly, according to the Proposition 13. II and the definition of (h n ,u n ), 
the following inequality holds 
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+ Khus^) + i' u " ■ Vc "ii^(-Bi, 1 ) + ii^rLi^) 

The estimate of F" is straightforward. From Lemma 12.61 we have 

ll^ll^(^) < C^^a^u^^ < CE\h\M\T). (4.9) 
With the help of Lemma 12.61 and interpolation arguments, we have 
K-Vcl^ ^||ul 4+S(4V) ||Vcl m 

^u^^Ju"^^ (4-10) 

In the same way, we can get 

II""' ' Vul^ Si) ^ CE?(h»,u n ,T). (4.11) 

From Lemma 12.61 Lemma 12.31 and interpolation arguments, we have 
- V/i n -Vu n 

h n Vh n ■ Vu" 
^||V/i n ■ V^Hii^) + 11—^7^—11^(5^) 

^c , l|v/i n || L oc ( ^ i) l|Vu"i| L ^ (i ji i) + ^11^ ^11^(4^)1^^11^^) ^ 

/ ft" 

^cII^IIl-^+^II^IIl^b^)^ 1 + II^^^IIl-^,) 
^c'll^Lo?^;^) KIL^) (i + II^ILs?^ 1 ,!), 

^CE 2 (h n , u n , T){1 + £(/i n , u n , T)). 



Similarly, we can get 
Vh n ■ D(u n ) 

M - hQ+ [ hn '- Wl^) < CE\h\M\T){l + E{h\u\T)). (4.13) 
Hence, from (14. 101) -( 14.131) . we gather 

ii u " • v^Wl^b^) + ii^rn^^j ^ 

^C{l + Afih 1 {l + E{h n ,u n ,T)))E 2 {h n ,u n ,T). 
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From (USD, (USD, (USD and ffl~T4"D . it follows 
IK*", ^ Ce A(52 ^[l + C7L 2 C< 2 (1 + W(l + A(7S(0)))S(0)]S(0). 
So we can choose £7(0) so small that 

A 2 



1 + CA 2 C 2 (1 + Afiho\l + ACE(0)))E(0) < 
e AC^(o) ^ ^±1 and 2C 1 ACE(0) ^ h Q , 



A + 2' 



(4.15) 



which yields || (h n , u n ) H^i ^ 4±±ACE(Q) for any T < f n . It follows that f n = 
T*. In fact, if f n < T*, we have seen that E{h n , u n , T n ) ^ ^fAC^O). So by 
continuity, for a sufficiently small constant a > we can obtain E(h n , u n , T n + 
a) ^ AC-E'(O). This yields a contradiction with the definition of T n . 

Now, if f n = T* < oo, then we have obtained E{h n ,u n ,T*) ^ ACE(0). 

As ||/i n || L f (B°< 1+e ) < 00 an d ll 11 "!!/, „CB°' e ) < °°j ^ implies that ||/i n ||L T ,(L2) < 
oo and Hu^U^^a) < oo. Thus, we may continue the solution beyond T* 
by the Cauchy-Lipschitz theorem. This contradicts the definition od T*. 
Therefore, the approximate solution (h n , u n ) ng ^ is global in time. 



4.3 Existence of a solution 

In this part, we shall show that, up to an extraction, the sequence (h n , u n ) ne pj 
converges in x M 2 ) to a solution (h,u) of ( 13.11) which has the desired 

regularity properties. The proof lies on compactness arguments. To start 
with, we show that the time first derivative of (h n , u n ) is uniformly bounded 
in appropriate spaces. This enables us to apply Ascoli's theorem and get 
the existence of a limit (h, u) for a subsequence. Now, the uniform bounds 
of the previous part provides us with additional regularity and convergence 
properties so that we may pass to the limit in the system. 

It is convenient to split (h n , u n ) into the solution of a linear system with 
initial data (h n , u n ) and the discrepancy to that solution. More precisely, we 
denote by (h^, u^) the solution to the linear system 

dthl + divu™ = 0, 

d t u n L - /iAu™ - 3/iV divu£ + f(u2) x + gVh n L - /3VAh n L = 0, (4.16) 

(^L> u i)i=0 = (h n ,U n ), 
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and (h n , u n ) = (h n - h% u n - u n L ). 

Obviously, the definition of (h n , u n ) entails 

h n —> ho — ho in B®'^ 6 , u n — > u in -B^'i' asn^ +00. 

The Proposition 13.11 insures us that 

(h n L , uj) -> u L ) in n (4.17) 

where u £ ) is the solution of the linear system 

d t h L + divu L = 0, 

d t u L - fiAu L - 3//Vdivu L + fuf + gVh L - [3V Ah L = 0, (4.18) 
(h L , u L ) t=0 = (/io - h , uo). 

Now, we have to prove the convergence of (^ n ,u n ). This is of course a 
trifle more difficult and requires compactness results. Let us first state the 
following lemma. 

Lemma 4.1 ((h n , u n )) neN is uniformly bounded in C^(M + ; fi 2 °i) x (^^ (^ + ! ^2 

Proof. Throughout the proof, we will note u.b. for uniformly bounded. 
We first prove that dth n is u.b. in L 2 (M + , 5° 1), which yields the desired 
result for h. Let us observe that h n verifies the following equation 

d t h n = -J n {h n divu") - J n (u" • Vh n ) - ho divu" + h divu£. 

According to the previous part, (/i n )„ g N is u.b. in L°°(B\^) and (u n ) ne pj is 
u.b. in L 2 (B\ l ) in view of interpolation arguments. Thus, — J n {h n divu") — 
J n (u n ■ V/i n ) — divu" is u.b. in L 2 (i?2 1 ). The definition of obviously 
provides us with uniform bounds for divuj in L 2 (i?2 1 ), so we can conclude 
that d t h n is u.b. in L 2 (B Q 21 ). 

Denote c n L = A" 1 divu£, c" = A" 1 divu", d n L = A" 1 div x u£ and d n = 
A -1 div^u". Let us prove now that <9 t c™ is u.b. in (L^r + L°°)(i?2 1 ) and 
that d t d n is u.b. in (Lr^ + L°°)(i?2 1) which give the required result for u n 
by using the relation u n = — A _1 Vc" — A~ 1 V ± d n . 

Let us recall that 

d t c n =AfiA(c n - c n L ) + f(<T - d n L ) + gA(h n - h n L ) + (3A 3 (h n - h n L ) 
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,._,,.,„_„ V/i n L>(u n )+V/i n divu" 
J n A 1 div u™ • Vu n - 2/i 



h n + h 

d t d n =u.A{d n - d n L ) - f(c" - c n L ) 

Vh n D(u n ) + V/i n divu 



- J n A~" div x u n • Vu n - 2/i 



h n + h 



Results of the previous part and an interpolation argument yield uniform 
bounds for u n in L 2 ^ (B £ 21 )nL^ (B^ £ ). Since h n is u.b. in L°° (B^), c n L and 
c n are u.b. in Lr^ [B\ 1 ), we easily verify that A(c n — c^) and J n A~ l div^u™ ■ 

Vu — 2/i v h^+h j are u m l-°2,iJ- Obviously, we have 

rf n+1 and u.b. in L 00 ^^). Because fr n and h n L are u.b. in Z 00 ^^), 
we have A(h n+1 - u.b. in L°°(i?2 -J. We also have h n and /i£ are 

u.b. in (Xi + L 2 ^ £ )(i?| 1 ) in view of Lemma 12.81 Thus, A 3 h n is u.b. in 
(Li +L 2 * £ )( J B° 1 ). So we finally get d t c n u.b. in (L 2 ^ + L°°)(5 2 ° 1 ). The case 
of dtd n goes along the same lines. As the terms corresponding to A 3 (h n — h?l) 
do not appear, we simply get d t d n u.b. in (L^ + L 00 )^^ -J. □ 

Now, we can turn to the proof of the existence of a solution and use 
Ascoli theorem to get strong convergence. We need to localize the spatial 
space because we have some results of compactness for the local Sobolev 
spaces. Let (Xp)peN be a sequence of C^°(1R 2 ) cut-off functions supported in 
the ball B(0,p + 1) of 1R 2 and equal to 1 in a neighborhood of -5(0, p). 

For any p G N, Lemma Hj] tells us that {{Xph n , Xp^ n ))n&i is uniformly 
equicontinuous in C(R + ; (-8° i) 1+2 ) an d bounded in E 1+£ . 

Let us observe that the application / i— > Xpf is compact from B^\ 
into B^i, and from 5f x into B 2 \. After we apply Ascoli's theorem to the 
family ((xph n , XpU n )) n£ N on the time interval [0,p], we use Cantor's diag- 
onal process. This finally provides us with a distribution (h, u) belong- 
ing to C(R + ;B 21 x (-B^'i) 2 ) and a subsequence (which we still denote by 
({h n , u n ) n£ ^) such that, for all p G N, we have 

(x P h n , X P u n ) -> (xph, as n -> +oo, 



in C([0,p]; J B° 1 x (B^) 2 ). This obviously infers that (/i n ,u n ) tends to (/i,u) 
in ^'( R+ x K 2 )- 

Coming back to the uniform estimates of the previous part, we moreover 
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get that (h, u) belongs to 
L°° (R+; B° 2 j +£ x (fijD 2 ) D L 1 (R + ; (Kg n £ 2 4 J e,3+£ ) x (flg* 6 ) 2 ) (4.19) 

and to C 1 / 2 ^; S° x ) x (C^(M+; B^)) 2 . 

Let us now prove that (h, u) := (/il,u^) + (/i, u) solves (13. ip . We first 
observe that, according to (14.11) . 

' h n t + J n (u n • Vh n ) + /i Ac" = - J n {h n divu™), 
u n + J n (u n • Vu") - /xAu" - 3/iV divu" + ^u") 1 + gV/i n - /3VA/i n 
W£>(u n ) + V/i n divu n 

The only problem is to pass to the limit in @'(R + x R 2 ) in the nonlinear 
terms. This can be done by using the convergence results stemming from the 
uniform estimates and the convergence results ( 14. 17ft and ( 14. 19ft . 

As it is just a matter of doing tedious verifications, we show, as an ex- 
ample, the case of the term V | ±™ ■ Denote L(z) = z/(z + ho). Let 
9 e C^(R + x M 2 ) and p G N be such that supptf C [0,p] x B(0,p). We 
consider the decomposition 

h 9Vh n divu" h 9Vh divu 
h + h n h + h 

=J n [6{l - L(h n )) Xp Vh n Xp div(u n L - u L ) + 9(1 - L(h n ))x P Vh n x P div(x P (u n ~ 
+ 9(1 - L(h n )) Xp V( Xp (h n - h)) divu + 9\7h Xp divu(L( Xp h) - L( Xp h n ))} 
, _ TN h 9Vh divu 

The last term tends to zero as n — > +oo due to the property of J n . As 
9L(h n ) and ft," are u.b. in L°°(i?2 i) and u£ tends to u^ in L X (B\ x ), the first 
term tends to in L X (B\ 1 ). According to (14. 1 91) . Xp( u " — u ) tends to zero in 
L 1 ([0,p]; -B 2 j) so that the second term tends to in L 1 ([0,p]; B\ 1 ). Clearly, 
X P ^ n Xpfr in ^(^.l) and L( Xp h n ) -> L(x p /i) in L 00 (L 00 D B^ x ), so that 
the third and the last terms also tend to in L^Z^i)- The other nonlinear 
terms can be treated in the same way. 

We still have to prove that h is continuous in B 2 \ +£ and that u belongs 
to C(M+; (Bl'lf). The continuity of u is straightforward. Indeed, u satisfies 

d t u = - u ■ Vu + fiAu + 3/iV divu - fu x - gV/i + /3VAft 
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h + h 

and the r.h.s. belongs to (L 1 + L°°)(i?2 i) by noting that we also have h G 

x ) in view of the interpolation argument. In a similar 
argument, one can obtain u G C(M + ; (5| l ) 2 ). We have already got that 
h G C(M+; Bl x ). Indeed, h -h e B^, u G L 2 (M+; , G 5^) 
and then o^/i G L 2 (M + ;i?2i) from the equation d t h = —h Q divu — div(/iu). 
Thus, there remains to prove the continuity of h in B\\ e . 

Let us apply the operator A k to the first equation of (13.1 j) to get 

A fc /i = - A fc (u • V/i) - ^ A fe divu - A fc (/i divu). (4.20) 

Obviously, for fixed the r.h.s belongs to Ll oc (R + ; L 2 ) so that each A^h is 
continuous in time with values in L 2 . 

Now, we apply an energy method to (I4.20p to obtain, with the help of 
Lemma 6.2 in fTUl . that 



lA^lli, ^ C||A fc / i |U 2 ( Q;fc 2- fc ( 1 + £ )||/ i ||^+e||u|| 42ii 



1 d 
2~dt 

+ ||A fc divu|| i2 + ||A fc (/idivu)|| L 2 
where Ylk a k ^ 1- Integrating in time and multiplying 2*( 1+ % we get 
2 fe ( 1+£ )||A^(t)|| L2 
^2 fc ( 1+£ ) || A fc (/i - ho) \W+C jf (a fc ||/i(T) bi+e ||u(r) 
+ 2 fc ( 2+£ )||A fe u(r)|| L2 +2 fc ( 1+£ )||A fc (/ i div U )(r)|| L2 )rfr. 

Since /i G L°°(^J e ), u G L\Bffi*) and /idivu G L^B^J 6 ), we can get 
]Tsup2 fc ( 1+£ )||A fc / i (t)|| L2 <\\ho ~ h \\ B iy + (l + \\h\\ LOC(B ^ ||u|| il( ^ ;2+e) 



+ 1 1 /i divu 1 1 x i < o°. 



Thus, Xjfcisgjv converges uniformly in L°°(R + ; B^) and we can conclude 
that h G C(R + ;Bl + 1 £ ). 
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4.4 Uniqueness 

Let (hi, Ui) and (h 2 , u 2 ) be solutions of (13. ip in E\ D with the same 
data (h — h , u ) constructed in the previous parts on the time interval [0, T}. 
Denote (5h, 5u) = (h 2 — hi, u 2 — Ui). From (13.11) . we can get 

d t Sh + u 2 • V5h + h div£u = F 2 , 

d t 8u + u 2 ■ V5u - fiAdu - 3/A7 div5u + f(5u) L + gV5h - /3VA5/i = G 2 , 

(5h,5u) = (0,0), 

(4.21) 

where 

F 2 = — 5u ■ V/ii — <5/i divu 2 — hi div<5u, 

„ r _ V<5/idivu 2 V/iidiv5u 

G 2 = - 5u • Vux + 2/x^= — + 2/i^= — 

ft + h 2 h + h 2 

+ 2 ^ ( TnTTT - ITTTT ) V/l1 divUl 
\h + h 2 h + hij 

+ VShD(u 2 ) + 2 V^Dgu) 

/l + ^2 ^0 + ^2 

\h + h 2 h + hij 
Similar to (13.11) . we can get 

\\(Sh,5u)\\ E , < Ce 011 " 2 "^'^^!!^!!^^^ + UGall^j). 

Noticing that h x G L^(B 2 j) n L^(B 2 j) and u 2 G L^(B% tl ), we can get 

11-^2 1 1 £1,(50 j <T= ll^uH^oo^o j 11x2(5^) + ll<^ILg?(B° J ll U 2 Hi^C-Bf^) 
+ IIMl~(B° J 11^1^(5^)- 



Moreover, from hi G L^(B 21 ) PI L^?(i? 21 ) by Lemma [2. 8 [ we have 

+ II^IL~(Bl jll^lll^BfJ + IIMz°?(B1 1 )II < Ml^(.B! 1 )- 
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Noting that hi, h 2 G L^(B^i), u 1; u 2 G L^B^i), and 

11^1 ||l°°([0,T]xR 2 ) ^ 2^°' II^2||l°°([0,T]xM 2 ) ^ 2^°' 

by the construction of solutions, we have 

11^211^(^0 J 

^II^IL-CBOjIIUllli^BfJ +4/i(l + ||'l2|| Z oo (A i i) )||<J/l|| £ o P( fli ii) ||U2|| i i, (A a i) 
+ 4/i(l + 11^211^(51^)11^111^(^^11^11^(^2,) 
+ 4 /ip/i|| L ^ (i ji a ) || /il || ^(Bi ^llui 11^(^2^) 

x (! + IML??^) + INL??^) + II^iIIl-^^II^IIl-^^))- 

Thus, we obtain 

||(5Mu)|k < Ce C||U2ll ^ (fl 2 2 1 »{ (l + T* + 4/^(1 + IIMx^;})) 

■II^iIIl ? > ( ^) + ^( t )}|I(^^)IL 

where limsup T ^ 0+ Z{T) = 0. 

Supposing that ^ 2 +^+ 2 (A + l)OE(O) < f besides f H45l) for E(0) and 
taking < T ^ 1 small enough such that C*|| u 2 1| z,i,(s| 1 ) ^ ln2 and Z{T) < |, 
we obtain \\(Sh, 5u)\\ E i^ = 0. Hence, (/ti,Ui) = (h 2 ,u 2 ) on [0, T]. 

Let T m (supposedly finite) be the largest time such that the two solutions 
coincide on [0, T m ). If we denote 

(hi(t), Ui(t)) := (hi(t + T m ), Ui(t + T m )), % = 1, 2, 

we can use the above arguments and the fact that 

H^ilU^OR+xR 2 ) ^ 2^° anC ^ II^*IIl°°(K+;_B§ jDB^ x ) ^ A(7.E(0) 

to prove that (hi,Ui) = (/i2,U2) on the interval [0,T m ] with the same T m as 
in the above. Therefore, we complete the proofs. 
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